High-fidelity fast quantum transport with imperfect controls 



M. Murphy/'Q L. Jiang, 2 N. Khaneja, 3 and T. Calarco 1 

Institut fur Quanteninformationsverarbeitung, Universitat Ulm, Albert- Einstein- Allee 11, 89081 Ulm, Germany 
2 Department of Physics, Harvard University, Cambridge, Massachusetts 02138 
3 School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138 

Effective transport of quantum information is an essential element of quantum computation. We 
consider the problem of transporting a quantum state by using a moving potential well, while 
maintaining the encoded quantum information. In particular, we look at a set of cases where the 
input control defining the position of the potential well is subject to different types of distortion, each 
of which is motivated by experimental considerations. We show that even under these conditions, 
we are able to perfectly transfer the quantum information non-adiabatically over any given distance. 
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In many current proposed implementations of quan- 
tum computers, it is necessary that we can transport 
computational states between operational sites. This 
state will often contain information from a previous op- 
eration that we want to preserve and use in a further 
operation. Consequently, we desire that the final quan- 
tum state becomes displaced but that the initial and final 
states are equivalent up to a global phase. In this sense, 
we will have preserved the encoded information. In prac- 
tice, this transport process is difficult to realise without 
altering the state: many external sources serve to heat 
or otherwise decrease the purity of the state while it is 
transported. In addition, as has been pointed out in 
transport processes may account for 95% of the operation 
time of a quantum computation. It is therefore advan- 
tageous to minimise the time required for this process 
while preserving the motional state. In this paper, we 
present an analytic solution for a one-dimensional sys- 
tem for transporting a quantum state over an arbitrary 
distance in the non-adiabatic regime using a harmonic 
potential, subject to a distortion of the input controls. 
This has application in systems where the control of the 
transport mechanism is imperfect, as is the case in many 
realistic experimental implementations, where a well de- 
signed input control becomes distorted, either through a 
limitation due to the experimental hardware, or through 
the interaction of the apparatus with an uncontrolled en- 
vironment. Analytic solutions to a driven quantum har- 
monic oscillator have been known for many years 0, Q , 
but there have been few attempts to utilise the results 
for the benefit of quantum information transfer Q . Our 
approach differs, since it is not focussed on any one im- 
plementation for transporting quantum states (although 
it has particular application to trapped particles), and 
we do not assume perfect control of the system. There 
are also some early attempts at high fidelity transport in 
experiments [H, H, 0, 0] . 

We first analytically solve the Schrodinger equation for 
a quantum state confined in a driven harmonic oscillator, 
and show that a suitable choice of the transport function 
d(t) results in an evolution of the wavepacket that con- 



stitutes what we shall refer to as 'perfect transport': the 
final evolved quantum state is equivalent to the initial 
state up to an irrelevant global phase which is analogous 
to a free evolution of the state in the frame of the poten- 
tial. This motivates us to treat the problem classically in 
order to derive a particular form of the driving function 
that satisfies this criterion. We examine the distortion of 
the ideal transport path analytically by the introduction 
of a general functional, and show that a large class of 
functionals that describe this broadening have no effect 
on the success of our transport. 

We model the transport of the particle from its ini- 
tial position over some arbitrary distance by a movement 
of the potential well according to the functional T>[d(t)] 
of an input function d(t) (our intended transport path), 
which defines the position of the centre of the potential 
well along the axis of transport for a given time t. The 
system can be modelled as a wave packet confined in a 
static harmonic potential of fixed frequency, subject to 
the Hamiltonian H(t) = \p 2 + \{x - V[d(t)}) 2 , with x 
and p being the usual quantum operators correspond- 
ing to position and momentum respectively, and where 
T>[d(t)] now plays the role of a driving function. Note 
that we have transformed the variables to make them di- 
mensionless. We prepare the system in a given eigenstate 
of the harmonic oscillator \ip n (x,t = 0)). If we denote 
the transport distance by Ax, our transport condition 
for the quantum case is that the fidelity between initial 
and final states be T = \("ip n (% — Ax, 0)\ip n (x, T))\ 2 = 1; 
in other words, the state be unchanged in the reference 
frame moving with the potential well except for a global 
kinematic phase. To verify this, we must solve the time 
dependent Schrodinger equation for our system. The nor- 
malised solutions are 0] 

t)) = e^ E » t+ ^- F '^f F \Mx, 0)) , (1) 

where E n is the energy eigenvalue corresponding to 
the nth eigenstate, and Tp is the translation operator 
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T F ip(x) = ip(x - F(t)). The function F(t) is defined as 



F(t) 



[ V[d(r)} sin(i - r)dr , 
./o 



and the phase <p(t) is given by 



<f>(t) = / 2 J F 1 (r) J F 1 "(r) + F" 2 (T)+ J F 1 ' 2 (r)dT 
Jo 
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We can satisfy our transport condition if the function 
F(t) satisfies F(T) = Ax, F(0) = F'(0) = F'(T) = 0, 
such that when t = 0, the wavefunction in equation (fT]) 
reduces to the normalised eigenstates of the harmonic 
oscillator, and at t = T, the wavefunction reduces to the 
original wavefunctions shifted by an amount Ax and with 
a global phase exp[— i(E n T + \(f>{T)}. 

There is a direct correspondence between the wave- 
function in equation {TJ) and the equations of motion 
for the classical analogue of our system, which can be 
seen by taking the expectation values of the operators 
x and p of the state in equation (TTJ), and then not- 
ing that the time evolution of these quantities obey 
the Newtonian equations of motion for a classical par- 
ticle confined within a harmonic potential well with con- 
stant frequency. The equations of motion are given by 
x(t) = p(i), p(t) = T)[d(t)] — x(t), where again all vari- 
ables are rescaled to make them dimensionlcss. Here, x(t) 
and p(t) refer respectively to the position and momentum 
of the classical particle along the axis of transport at time 
t. These equations have solutions x(t) = x c (t) + F(t) 
and p(t) = p c (t) + F'(t), with F(t) as defined in equa- 
tion ([2]). The functions x c (t) and p c (t) are the solutions 
to the homogeneous equations of motion, which therefore 
describe the simple harmonic motion undergone by the 
particle when no transport is undertaken. In the classical 
picture, the condition for performing perfect transport 
becomes x(T) = x c (T) + Ax and p(T) = p c {T), where 
Ax denotes the displacement of the potential well at the 
final time. Our perfect transport condition specifies the 
boundary conditions on F(t) and its first-order deriva- 
tives. We can rewrite F(t) in terms of T>[d(t)] by noting 
that equation ([2} is a Voltcrra integral equation of the 
first kind with a trigonometric kernel Q, and hence has 
the solution 



V[d(t)]=F"(t)+F(t) 



(4) 



We first search for a form for the function d(t) by sup- 
posing that V[d(t)] = d{t). Substituting the transport 
conditions on x(t) and p(t) into equation (J5]) and impos- 
ing the condition F"(0) = F"(T) = 0, we find F(0) = 
and F(T) = Ax. If we fix the condition that d'(t) = 
when t = and T, we subsequently place boundary con- 
ditions on F'(t) and F"'(t) at the initial and final times. 
Taking all of these conditions into account for F(t) and 
its derivatives, we can now construct a function d(t) that 



transports our particle perfectly. This procedure is as fol- 
lows: we construct a general F (t) by taking the simplest 
form of transport function (a linear function) and adding 
a series of Fourier components. We scale the components 
so that their period matches the transport time. We then 
apply the boundary conditions on F(t) to solve for the 
Fourier coefficients. Due to the periodicity of the com- 
ponents, we have only five independent boundary con- 
ditions, and so we may uniquely specify only this many 
Fourier components. Substituting this into equation ((4|) 
gives us the solution 
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which we shall call our transport function. This equa- 
tion depends on the transport time T, so that one must 
choose the correct transport function for the appropriate 
transport time. Since we have scaled all the variables in 
our system, T = lix represents transport over one pe- 
riod of the harmonic motion in the potential well (the 
trap period). We henceforth assume that d{t) = d a (t). 
(It should be noted that one can follow a more rigorous 
derivation by starting from a description in the frame of 
an optimal control problem, but this was not presented 
here. See for instance 0.) 

Through consideration of the above, we may state 
the following: if the functional T>[d(t)] has a form such 
that the boundary conditions on F(t) are preserved (and 
V[d(t)] is non-singular for all t), then the functional V 
does not affect the transport of the particle. We consider 
the following three forms for T>[d(t)], and briefly discuss 
the motives for doing so. 

(i) The d(t) model: V[d{t)} = d(t) + ad(t), where a is 
a real constant, and d(t) represents differentiation with 
respect to time. A physical interpretation for this model 
could be that we consistently 'overshoot' or 'undershoot' 
our desired potential well position, so that as we move the 
well more quickly, the deviation from the desired position 
becomes greater. 

(ii) The piecewise model: The functional T> casts d(t) 
into a piecewise form 



V[d(t)} = d(t n ) for te[t n -^,t n + ^ 



(6) 



where t n = ^ for a given JV e Z+, N > 1. This has 
the effect that the potential undergoes discrete 'jumps' 
in its position along the transport axis, which could be 
due to a sampling rate limitation of the equipment used 
for experiment. 

(iii) The Fourier model: T>[d{t)\ = d(t) + g(t), where 
g(t) is a discrete Fourier series. This is of relevance since 
we may decompose any periodic signal (for instance, pe- 
riodic pulse distortions) into a Fourier series (note that 
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FIG. 1: (Colour online) Distortion of the transport function 
over time period T = 2n. The thick (red) solid line is the 
transport function from equation (0, dot-dashed (green) line 
is the d(t) model (a = 1), the dashed (blue) line is the Fourier 
model (random coefficients in the range [-1,1]), and the dotted 
(black) line is the piecewise model with an exponential-type 
smoothing (N — 8). 



the Fourier representation must converge at all times t 
to the signal being represented). 

We now show that each of these functionals the condi- 
tions for perfect transport are still satisfied. In the case of 
the d(t) model, we begin by writing d(t) = F"(t) +F(t). 
Hence V[d(t)] = F"(t)+F(t), where F(t) = F'(t)+F(t). 
It should be clear that F(t) satisfies the same boundary 
conditions as F(t). Hence we conclude that we achieve 
our desired transport with this functional for any trans- 
port time T. 

In the case of the piecewise function, we substitute the 
functional directly into equation and solve for F(t). 
We find in the limit that T -» T k = 2kir, fce2+ the 
boundary conditions for F(t) are satisfied as before. In 
other words, the piecewise functional can only achieves 
the boundary conditions when the transport time is an 
integer multiple of the period of the harmonic trap. It 
may additionally be the case that the movement of the 
potential is not exactly stepwise, but that instead the 
movement is smoothed out (for instance, if we consider 
a segmented ion trap, this will be due to the charging 
characteristics of the electrodes). We can model this by 
writing equation ((6]) as 

V[d(t)] = d(t n ) - q(t - t n + ^) [d(t n ) - d(i„_i)] (7) 

so that q(t) describes the smoothing from the previous 
value in the stepwise function to the next. Substituting 
this into equation j2]), we can calculate that the part of 
the integral dependent on q(t) evaluates to zero. Hence 
we may conclude that any smoothing of the transport 
path due to these terms may be ignored. 

Finally, we consider the Fourier model. Again, through 
direct integration of the functional via equation ([2|), we 
can obtain the associated function F(i). Here, we see 
that we satisfy the boundary conditions if the period of 
the function g(t) is Tfc/2, with Tk as given above. We can 
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FIG. 2: (Colour online) Trajectories through phase space in 
the frame of the moving potential of a transported particle 
over time period T = 2n. See caption of Figure fT] for key. 
For illustative purposes, the particle was given some initial 
momentum. The centre circle is the trajectory of a freely 
oscillating particle with constant energy and the same mo- 
mentum in a static potential. 
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FIG. 3: (Colour online) The evolution of the ground state 
probability distribution from (a) an initial time t = to (f) 
the final time T — 2tt. See caption of Figure [T] for key. 



satisfy this by tuning the the frequency of the harmonic 
potential to accommodate the periodic noise. Figure [1] 
shows sample transport paths 2?[cf(t)] over a short trans- 
port time T = 2tt. One can see that the deviation here is 
not small; we significantly disturb the motion of the par- 
ticle. Figure [2] shows the classical trajectories through 
phase space of the particles transported according to the 
transport paths T>[d(t)] (or, equivalently, the expectation 
values of the operators x and p). One sees here that the 
trajectories begin and end on the constant energy curve 
given by the free oscillation in the well. 

Figure[3]shows snapshots of the evolution of the ground 
state wavefunction subject to distortion; although the 
probability distributions diverge from each other at inter- 
mediate times for the different models for T>[d(t)], at the 
final time the wavefunctions converge to the initial pro- 
file of the ground state probability distribution displaces 
to the final position. Furthermore, one may verify that 
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FIG. 4: (Colour online) The fidelity between the instanta- 
neous ground state and the transported ground state over a 
time T — 2n. See caption of Figure [1] for key. 

the energy expectation at the final time is the predicted 
value of (E) = E n = (n + \) 1 which can be immediately 
found from equation JT]) (although (E) is only uniquely 
specified for the ground state). 

The scheme also allows us to transport far from the 
adiabatic limit, which can be demonstrated by calculat- 
ing the fidelity between the actual transported ground 
state and the instantaneous ground state of the displaced 
potential well for different values of the transport time. 
As we transport the state, the fidelity deviates further 
from unity at intermediate times, only to recover again 
at the final time (demonstrated in Figure I3J . Although 
here we have chosen the ground state of the potential, 
we may well have used any of the other eigenstates of 
the harmonic oscillator and produced similar results. 

We can also begin with superpositions of states. Con- 
sider the coherent superposition of eigenstates given by 
^(XjO)) = J2 n =o CnlV'n^i 0)), where c n are normalised 
coefficients. In order to recover the maximum fidelity, 
we must preserve the relative phases between initial 
states. The nth eigenstate acquires a relative phase 
exp[— i(n — m)T] with the mth eigenstate during trans- 
port, and so in order to preserve its phase relation with 
the other superposed states, we can choose T = 
Hence the transport condition is satisfied. 

If we instead begin with a mixed state described by 
the density matrix p{x, 0) = X)«Lo PnlVAi) (VVi I, and take 
the fidelity^ T = (tr[^pTp 2 *Jpi\ l/2 ) 2 where pi and 
P2 are the initial and final states respectively, we can 
show that the fidelity of transport does not depend on the 
coefficients p n , so that the distribution of states remains 
constant during transport. We can then infer that the 
transport is insensitive to temperature. 

In this paper, we have derived analytic solutions for 
the transport of a quantum state via a moving harmonic 
potential with a constant frequency. We considered the 
conditions under which a functional that distorts the in- 
put control achieves the conditions we have set for per- 
fect transport. In particular, we studied three different 
models that have a quantitative relevance in experiments 
dealing with quantum transport, and showed that under 



certain conditions, all three models describe a broadening 
of the transport path without detriment to the transport 
success. We briefly review these conditions, (i) The d(t) 
model fulfills the transport condition for any transport 
time, (ii) The piecewise model fulfills the transport con- 
dition when the transport time is an integer multiple of 
the trap period, (iii) The Fourier model (which models 
signal distortions) satisfies the boundary conditions if the 
signal is periodic with half the trap period. 

Of course, in a realistic situation beyond our harmonic 
oscillator model, complete insensitivity to such a broad 
range of control imperfections is not to be expected: in 
that case, we not only have distortion of the input con- 
trols, but also distortion of the shape of the potential 
itself. However, our result indicates significant robust- 
ness (by which we mean a low sensitivity of the trans- 
port fidelity to such distortions) is likely to be obtained 
at least when potential anharmonicities are small, since 
the system is well approximated by the harmonic oscil- 
lator. This will be the subject of future investigations. 
Such deviations may also be overcome by application of 
optimal control methods [HI, [12]. The consequence of 
performing robust transport, particularly in the presence 
of such imperfections, is that we can distribute quantum 
information to separate elements of a quantum computer 
with a very high fidelity in short times and over large dis- 
tances. This is essential for scalability and fault tolerance 
of quantum systems for future use in computation. 
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